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Abstract 

^ ' The chiral-odd structure function and the tensor charge of the deuteron are 

' studied within the Bethe-Salpeter formahsm for the deuteron ampHtude. Utihzing a 

' simple model for the nucleon structure function, hi ^ hi{x) is calculated and the nuclear 

, effects are analyzed. 

0^ 



"o. 1 Introduction 



Oh! 

1^ ■ There is an increasing interest in the transverse quark distributions in the nucleon, especially the 

chiral-odd structure function of twist two [|T], 0|, and its first moment, the tensor charge 

\ of the nucleon p, ^, ||, H, 0, |^. The reasons for the interest are the following. First, hf{x), 

^ ! complimentary to the conventional structure functions ^'^d gi2j is essential to understand 

the nucleon spin substructure. At the same time, the tensor charge, g^ , like the axial charge, 
is one of the fundamental observables for the nucleon. So far, however, little is known about 
both hi{x) and g^ . Second, semi-inclusive deep-inelastic experiments aiming to measure , 
and hence the tensor charge of the proton and neutron, g^ , are being planned |^. Now is the 
time for the theory to make predictions. 

It is important to remember that experimental results for the neutron can only be obtained 
from the indirect experiments, using the the lightest nuclear targets, the deuteron and ^He. 
Therefore, it is necessary to have a reliable estimate of the nuclear effects, in order to extract 
dependable information about the neutron from the nuclear data. 

^INFN Postdoctoral Fellow 



1 



In this letter we study the chiral-odd structure function and the tensor charge of the 
deuteron, h^{x) and g^, in the Bethe-Salpeter formahsm. The main motivation is to study the 
nuclear effects in the transverse-polarized deuteron structure function, /if, and the deuteron 
tensor charge, g!^. Since the nucleon tensor charge measures the net number of transversely 
polarized valence quarks (quarks minus antiquarks) in a transversely polarized nucleon , and 
so it gets no contribution from the sea-quarks, the neutron isoscalar tensor charge is equal to 
the one of the proton. As a result, the measurement of the deuteron tensor charge may, in prin- 
ciple, provide a test for the models of the nuclear effects if the proton tensor charge is known 
precisely. In addition, an independent measurement of the neutron's h"^ and g^ is obviously 
useful to understand the transverse spin distribution in the nucleon. Therefore it is essential 
to account for the nuclear effects in the deuteron for both the chiral-odd structure function /if 
and the tensor charge, g^. 



2 Basic formulae 

The chiral-odd structure function of the deuteron of twist two, /if (x), is defined as a simple 
extension of the case of spin-half hadrons . It can be expressed in terms of a matrix element 
of the quark bilocal operator: 
f fix 

/ — e^^^(PD Sd I ^(0)i(T'^^75^(An) | Pd Sd) = 2/if (x) ■ T"" + (higher twist terms). (1) 
An antisymmetric tensor, T^'^, in the rest frame of the deuteron is defined by: 

T^^ ^ (S^P^ - S^P^) , (2) 



where P^ is the deuteron momentum, Sj^ is the spin of the deuteron 



SUM) = -^e^-^^E:iM)Ep{M)Pn„ (3) 

where E{M) is the deuteron polarization, with M = ±1, 0, Mo is the deuteron mass. So 
satisfies the conditions. Si) = —1 and PdSd = 0. 

The tensor charge of the deuteron is the first moment of /if (x) and can be readily obtained 
from eq. (|l|): 

{Pd Sd \ mta'^'^^.m \ Pd Sd) = 2g^{x)T^\ (4) 
Using the following identities in the rest frame of the deuteron: 

T'''T,,\m=i = -2M|, (5) 

zT^V^,75Im=i = -2Mb75737o, (6) 
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the structure function /if (x) and the charge are rewritten as the matrix elements of the 
Lorentz scalar operators: 

h^x) = ^J^e^'^Pr,SD\miBl3lo^{Xn)\Pj,So), (7) 

9? = -^{PDSD\miBl3lom\PDSn), (8) 

It is impossible to solve eqs. (|^) and (||) directly, since the hadron state, | Pd Sd), still cannot 
be described in QCD. Instead, we use approach where the "elementary" structure function, 

, is obtained on the quark level, while the nuclear effects are calculated on the hadron 
level (e.g. [jTT], |T^] and references therein). We utilize the Bethe-Salpeter formalism, which was 
previously applied for studying the structure functions F^2: ^1,2 (^f [JTT], |T^. In what follows, 
we neglect small possible off-mass- shell effects in the nucleon structure functions [0, P^ . 

In terms of the Bethe-Salpeter amplitudes /if reads: 

,D( [ d^P .n( a;m \ Tr{^fA/(po,p)75737o^A/(po,p)(P2-m)}|^^^^ 

[X) — I I hi ; — - ; r , (9) 



(27r)4 \Pw + Pi3j 2(pio-hpi3) 
where m is the nucleon mass, p is the relative momentum of nucleons in the deuteron, piQ 
and pi3 are the time and 3-rd components of the nucleon momentum, pi, and pi^2 = -Pd/2 ±p. 
^Af(PO) p) is the Bethe-Salpeter amplitude for the deuteron with the spin projection M and the 
isoscalar structure function /if is defined as /if (a;) = {hi{x) + /i"(x))/2. Note, that there is no 
explicit mesonic (sea quarks) contribution to the structure function /if. The mesons manifest 
themselves by binding the nucleons in the deuteron and, therefore, defining the structure of the 
amplitude \l'jvf (po, p). 

Then eq. (^) can be rewritten in the convolution form p7[ ]: 

Mo/m 

/if(x)= / |/if (^^j /7^/^(y), (10) 

where H{y)^/^ is the "effective transverse distribution" of nucleons in the deuteron and its 
definition is obvious from eq. (^. This is a usual result for the nuclear structure functions 
calculated in the leading twist approximation, neglecting the off-mass-shell deformation of the 
nucleon structure function. Eqs. @) and ( p!o|) allow immediately to write the sum rule for the 
deuteron tensor charge: 



9? 



Molm 

dxh^ix)- f dyH^'l^iy) (11) 



g'^ ■ {Pd Sd \ iV(0)75737oiV(0) | Pd Sd), (12) 
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where N{x) is the nucleon field. Eqs. (|TT]) and (|I2|) define the renormahzation factor of the 
nucleon tensor charge by the deuteron structure. This sum rule is similar to the sum rule for 



the deuteron structure function, [13 



± 1 

J dxng^ixD) = J dxg^ix) ■ {Pd Sd \ N{0)'y,^sN{0) \ Pd Sd). (13) 







This sum rule is used for estimating the integral of the neutron structure function, g^, from 
the combined proton and deuteron data. 

Note that in the static limit the matrix elements on the r.h.s. of eqs. (0) and (p!3|) are equal 
and approximate formulae for the deuteron structure functions are: 

gn^D)^{l-lwn)gnx), (14) 
h?{xn) ^ (l - Iwn) <(x), (15) 
where wd is weight of the d-wave component in the deuteron wave function. 



3 Numerical calculations 



The method to calculate numerically expressions like (^ is discussed in ref. |T3[. The most 
important details of the calculations are: 

1. A realistic model for the Bethe-Salpeter amplitudes is essential for a realistic estimate of 



the nuclear effects. We use a recent numerical solution [11| of the ladder Bethe-Salpeter 



equation with a realistic exchange kernel 16 



2. The Bethe-Salpeter amplitudes and, therefore, eq. (^ have a nontrivial singular structure. 



These singularities must be carefully taken into account |18, 19, 20 



3. The BS amplitudes are numerically calculated with the help of the Wick rotation. There- 
fore, the procedure of the numerical inverse Wick rotation must be applied. 

4. An essential ingredient of the calculations, the nucleon structure function hf{x), is un- 
known both theoretically and experimentally. A reasonable estimate for using in our 
calculations should be found. 

The effects of the nuclear structure, the Fermi motion and binding of nucleons, is encoded 
in the effective distribution function H^^^{y) in eq. (|10D. The result of calculation of H^''^{y) 
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with the Bethe-Salpeter amphtude of the deuteron ||Tl], |1^] is shown in Fig. (sohd hne). Two 
other effective distributions for the deuteron are also shown in Fig. ^ for comparison. The 
first (dashed hne) is f'^^^iy), which defines the spin-dependent structure function g^. Exphcit 
expression for f^^^{y) is similar to the one for Hi^^{y) (see eqs. (^ and (0)), but with 
the operator 75(70 + 73) instead of 757370- If we again replace the operator by (70 + 73), we 
get another distribution, f^^^{y) (dotted hne), which defines the spin-independent structure 
function F2 . Each of the three effective distribution functions {Hi^^{y), f^^^{y) or f^^^{y)) 
would represent the correspondent structure functions of the deuteron (/if, or F^) if the 
nucleons in the deuteron would be the elementary fermions. 

Note that the function H^^^{y) is smaller at the maximum than f'^^^{y), and also slightly 
smaller in its normalisation: 



J dyH'"''{y) = {Pd Sd \ iV(0)75737oiV(0) | Pd Sd) = 0.9208, (16) 


J dyf^'iy) = {Pd Sd \ N{0)^,jsN{0) \ Pd Sd) = 0.9215. (17) 



This is opposite to the case of the nucleon, where hi is expected to be larger than g^ (after 
exclusion of the quark charge factors). The effect is different, since in the polarized deuteron 
both nucleons are essentially polarized along the same direction, while in the polarized nucleon 
one quark is polarized in opposite direction and somehow cancels contribution of one of the 
quarks polarized along the nucleon polarization. The difference between H^^^{y) and f^^^{y), 
and between the tensor and axial charges, (0) and (|17D, is extremely small since the deuteron 
in essentially a nonrelativistic system and matrix elements of the operators 7573 and 757370 are 
exactly the same in the static limit. 

The numerical value in eq. (|l^) is the deuteron structure factor renormalizing the nucleon 



tensor charge in the deuteron, calculated in a particular model of the A^A^- interaction ||T6 | 



Taking into account the differences among existing models, we can estimate a possible "model" 
error in (|T§): {Pd Sd \ iV(0)75737oA^(0) | Pd Sd) = 0.93 ± 0.015. 

To calculate the realistic structure function h^lx) we need the nucleon structure functions 
hi{x). However, so far there is no existing experimental data for this function, and very little 
is known about the form of in theory. In the present paper we follow the ideas of ref. to 
estimate . Since the sea quarks do not contribute to , its flavor content is simple: 

h^{x) = 6u{x) + 6d{x), (18) 
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where 6u{x) and 6d{x) are the contributions of the u- and d-quarks, respectively |^ ^. Since 
the matrix elements of the operators oc 7573 and oc 757370 coincide in the static limit, as a 
crude estimate we can expect that 

Su{x) ~ Au{x), Sd{x) ~ Ad{x), (19) 

where Au{x) and Ad{x) are contributions of the u- and d-quarks to the spin of the nucleon, 
which is measured through the structure function g'^^ . Subsequently, the simplest estimate for 

h^{x) = aAu{x) + (3Ad{x), (20) 
with a = P = l (21) 

should not be too unrealistic. In fact, the bag model calculation shows that difference between 
6q and Aq is typically only a few percent [||. This analysis is mostly a qualitative one, since it 
is limited by the case with one quark flavor and does not pretend to describe phenomenology. 
To evaluate possible deviations from the simplistic choice of , (pO|) with (pi]), we suggest: 



a = 5u/Au, P = Sd/Ad, (22) 

where 6q and Aq are the first moments of 6q{x) and Aq{x), respectively (g = u,d). For 6u and 
6d we can adopt the results from the QCD sum rules and the bag model calculations P, ^. As 
to Am and Ad, we can use the experimental data analysis or theoretical results, e.g. the 
QCD sum rules results p2[. Thus, we estimate 



a = 1.5 ±0.5, /? = 0.5±0.5, (23) 

at the scale of = 1 GeV^. 

The realistic form of the distributions Au{x) and Ad{x) can be taken from a fit to the 
experimental data for . In our calculations we use parametrization from ref. |]23[. At this 
point we have to realize that, in spite of the expected relations (|T^), distributions 6q and Aq 
are very different in their nature. Especially at x < 0.1, where Aq probably contains a singular 
contribution of the polarized sea quarks, but 6q does not. Therefore we expect eq. ( ppj ) to be 
a reasonable estimate in the region of the valence quarks dominance, say a; > 0.1 (see also ||^). 
For completely consistent analysis, the parameters a and (3, and the distributions Au{x) and 
Ad{x) should be scaled to the same value of Q^. However, for the sake of the unsophisticated 
estimates we do not go into such details. 
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The results of calculation of the nucleon and deuteron structure functions, (solid lines) 
and /if (dashed lines), are shown in Fig. 0. The group of curves 1 represents case (pT]), which 
is a possible lower limit for hi'^ in accordance with our estimates (|23D. Curves 2 represent the 
case a = 1.5, (3 = 0.5, which is close to the mid point results of the bag model and the QCD 
sum rules. The upper limit corresponding to the estimates (p3D is presented by curves 3. For 
all cases the deuteron structure function is suppressed compared to the one for the nucleon, 
mainly due to the depolarization effect of the D-wave in the deuteron. This is quite similar to 
the case of the structure functions and g^. To illustrate this similarity we show the ratio 
of the structure functions h^/h^ (solid line) in Fig. ^ together with the ratio gi / gi (dotted 
line). The straight dash-dotted line in Fig. ^ represents the value of the matrix elements ([16D 
and (|17D , which are not distinguishable on this scale. This line approximately corresponds to 
eqs. (|1D and (||). 

Note that our estimate of the nucleon structure function h'^ , (^), gives systematically a 
larger value for the function than naive suggestion (^), the curves 1 in Fig. ^ which essentially 
corresponds to the estimate hi ^ (18/5)5'^, neglecting possible negative contribution of the s- 
quark sea [^. The large size of the effect suggests that it can be detected in future experiments 
with the deuterons 0. The ratio of the structure functions /h^ in this case is also shown 
in Fig. ^ (dashed line). 

4 Summary 

We have considered the chiral-odd structure function of twist two, /if (x), and its first moment, 
the tensor charge of the deuteron. In particular, 

1. The structure function of the deuteron is calculated in the Bethe-Salpeter formalism 
for the deuteron amplitude. Explicit analytical and numerical results for the effective 
distribution function of the nucleons in the deuteron, defining the deuteron structure 
function /if are obtained. 

2. It is shown that the deuteron structure function /if (x) exhibits deviation from the nucleon 
structure function /if similar to the case of the usual spin-dependent structure functions 
g^ and g^ . The deuteron structure factor renormalizing the tensor charge of the nucleon 
in the deuteron is calculated. 

3. A simple phenomenologically motivated model for the evaluation of the nucleon function 
/if is suggested, using the available results of the QCD sum rules and the bag model for 
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. The estimates show that the nucleon structure function, hi{x), can be significantly 
larger than the spin-dependent structure function after extraction of the averaged quark 
charge, i.e. {18/5)gi{x). This effect definitely can be tested experimentally in the semi- 
inclusive deep inelstic scattering and, therefore, can bring a better understanding of the 
spin structure of the nucleon. 
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Figure captions 



Figure 1: The effective distribution function H^^^iy) (solid line) for the nucleon contribution to 
the deuteron structure function /if (x). Other effective distributions: spin-dependent f^^^{y) 
(dashed hne) and spin- independent f^^^{y) (dotted line). 

Figure 2: The structure functions h^{x) (solid lines) and h^{x). Groups of curves: 1 - a = 
1.0, P = 1.0 ("lower limit"); 2 - a = 1.5,/? = 0.5 ("mid point"); 3 - a = 2.0,/? = 0.0("upper 
limit") (see eqs. and (|23|)). 

Figure 3: The ratio of the deuteron and nucleon structure functions h^{x)/h^{x). Curves: 
solid line - a = 1.5, [3 = 0.5; dashed line - a = 1.0, (3 = 1.0. The dotted curve presents the ratio 
of the spin-dependent structure functions gi{x)/g^{x). The dash-dotted line approximately 
corresponds to eqs. (0) and (|15|). 
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